Non-resonant inelastic x-ray scattering involving excitonic excitations 
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In a recent publication Larson et al. [ij] reported remarkably clear d-d excitations for NiO and CoO 
measured with x-ray energies well below the transition metal K edge. In this letter we demonstrate 
that we can obtain an accurate quantitative description based on a local many body approach. We 
find that the magnitude of q can be tuned for maximum sensitivity for dipole, quadrupole, etc. 
excitations. We also find that the direction of q with respect to the crystal axes can be used as an 
equivalent to polarization similar to electron energy loss spectroscopy, allowing for a determination 
of the local symmetry of the initial and final state based on selection rules. This method is more 
generally applicable and combined with the high resolution available, could be a powerful tool for 
the study of local distortions and symmetries in transition metal compounds including also buried 
interfaces. 

PACS numbers: 78.70.Ck, 78.20.Bh, 71.70.Ch 
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Transition metal compounds with partially filled d- 
shells show a large variety of interesting properties, like 
metal insulator transitions.colossal magneto resistance 
and super conductivity 0, y| • One reason for the com- 
plex behavior of these materials is the strong interplay 
between the orbital, charge, spin, and lattice degrees of 
freedom In the manganates, for example, the orbital 
and charge ordering could be related to the magneto re- 
sistance [5| . For vanadium oxides it has been shown that 
the orbital occupation changes drastically at the metal 
insulator transition 0, 0] and in the layered high tem- 
perature superconducting cuprates the two dimensional 
electronic structure is intimately linked to the unoccu- 
pied dx-2_yi orbital. More recently there is a lot of activ- 
ity in using interface induced effects in transition metal 
oxides to strongly modify the physical properties Q . 

Since the orbital degrees of freedom play an impor- 
tant role in all of these materials it is highly desirable 
to have good experimental methods to determine the en- 
ergy scale of the crystal or ligand field splitting and the 
local symmetry. In principle this can be done by optical 
spectroscopy [lO| , however these so called d-d excitations 
are even and therefore optically forbidden and often com- 
pletely masked by transitions involving small amounts of 
impurities or defects. The reason why some of these tran- 
sitions are optically visible at all is due to simultaneous 
excitations of magnons or phonons. This results in an 
intensity typically 1000 times smaller than the intensity 
found for the close-by charge-transfer or Mott-Hubbard 
excitations. In multilayers and interfaces the problem is 
even more severe since there may be a variety of optical 
transitions due to other components which quickly mask 
out the d-d transitions. One is even not always able to 
easily discriminate between absorption peaks due to d- 
d excitations also referred to as orbiton excitations and 
multiple phonon excitations for example 



techniques have been developed to study d-d excitations. 
At the transition metal if -edge a Is to 4p excitation is 
involved. In the intermediate state the d levels shift due 
to the changed local potential. This energy shift can 
change the occupied d orbital wave function, which may 
result in a resonant enhancement of the d-d excitations 
11, Q El • The draw-back of resonant scattering at the 
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Recently resonant inelastic x-ray scattering (RIXS) 



i^T-edge is however, that charge-transfer excitations are 
enhanced much more efficiently than the d-d H, 3] , for 
two reasons. The first is that the 4p-orbitals of the in- 
termediate state are quite spatially extended and have a 
small interaction with the 3c? orbitals but a very large one 
with the surrounding O 2p orbitals. The second reason 
is that the spherical core hole potential does not enhance 
d-d transitions directly. An other option developed re- 
cently is RIXS at the transition metal L2.3 edge or M2,^ 
edge. Here one excites (and de-excites) a 2p or ip tran- 
sition metal core electron into the 3d valence shell. With 
this technique one can choose which of the low-lying en- 
ergy states one wants to enhance by selecting the incident 
energy and polarization ll, 

A recent paper by Larson et al. l[ exhibited clear d- 
d excitations within the gap in NiO and CoO with the 
use of non-resonant inelastic x-ray scattering (NIXS) for 
energies just below the K edge. In principle one should 
expect that these excitations can be seen with NIXS, but 
surprisingly they found that the intensity of the d-d ex- 
citations at certain q vectors is much higher than the 
intensity of the charge-transfer or Mott-Hubbard excita- 
tions. They analyzed their experimental findings in the 
framework of LDA-I-U, which describes the transitions in 
terms of a one particle interband transition rather than 
a transition involving strongly bound excitonic states as 
is known to be the case for these states in NiO and CoO. 
Very interesting to note though is that the angular de- 
pendent results can be quite well described within the 
band structure approach for the case of NiO because the 
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transition involves basically a promotion of a t2g elec- 
tron into an unoccupied eg state. In more complicated 
cases involving multi Slater determinant excitonic bound 
states the situation will be quite a bit more involved as 
also realized by the authors of that paper and previous 
work on the cuprates by Ku et al. [20| • 

In this letter we will develop a local but many body 
treatment of NIXS and describe the observed d-d exci- 
tations within a configuration interaction cluster calcu- 
lation analogous to the approaches used for analyzing 
the energy positions of optical d-d excitations [2l[ or the 
RIXS spectra at the L2,3 edge [H, The goal is to 
show that we can have a straightforward and quantitative 
description of the NIXS process and that this will open 
up new opportunities to extract detailed and invaluable 
information concerning the local electronic structure of 
correlated electron systems not easily accessible by other 
techniques. 

The interaction of matter with light is given by two 
terms. One proportional to the vector potential {A) 
squared, the other proportional to the dot product of 
the momentum operator for the electrons (p) with the 
vector potential. 



Hi, 



A' 



-p • A 



(1) 



At resonance the second term {p ■ A) is responsible for the 
largest contribution to the scattering cross section. Off 
resonance however this term looses importance rapidly 
and the scattering is governed mainly by the term A^. 
The off resonance scattering cross section is then given 
by 
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One can define the dynamical structure factor S{q^ lj) 
which is a function of the scattering vector q — ki - 
kf and the energy loss uj = tOi — ujf as S{q,uj) = 



-./{rl 



), which has the advantage that all 
non-material dependent factors arc factored out. 

S{q, u)) is a sum over transition probabilities multiplied 
by a delta function responsible for the energy conserva- 
tion. This can be written as a Greens function in the 
spectral representation: 

S{q,Oj) = Y.VfW^'^ 
f 
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With the transition matrix equal to T = e* ' . 

To enable a rather direct symmetry analysis we prefer 
to discuss the transitions in terms of monopole, dipolc, 




FIG. 1: (color online) Expectation value of the spherical 
Bessel function for monopole, quadrupole and hexadecapole 
transitions as a function of q. 



quadrupole, etc. excitations and in order to do so we 
expand the transition matrix on spherical harmonics. 
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with Cm' = 2h+i ^knn and Yk.m thc spherical harmon- 
ics. This results in a sum over fc of a spherical Bessel 
function of order k times a spherical harmonic of order 
k. For d-d excitations only monopole (fc=0), quadrupole 
{k=2) and hexadecimalpole (A:=4) transitions are allowed 
and therefore only three values of k have to be evaluated. 
Wc first discuss the effect of changing the length of the 
q vector, which enters via the expectation value of the 
spherical Bessel function over q r and second discuss the 
effect of changing the orientation of the sample with re- 
spect to the q vector which enters as a spherical harmonic 
over the angular coordinates of q. 

The expectation value of a spherical Bessel function 
becomes small if the spherical Bessel function oscillates 
many times over the length scale of the product of the 
initial and final state wave function. We therefore expect 
the maximum intensity to occur at q values correspond- 
ing to thc period of a spherical Bessel function of length 
comparable to the atomic radial extent of the d wave 
function. In figure 1 we plot the expectation value of 
thc spherical Bessel function for k=0, 2 and 4 calculated 
for Ni^+ where the radial wave function has been cal- 
culated within the Hartree-Fock approximation with the 
use of Cowan's code (2^. One can clearly see that each 
of the different multipoles has a maximum at a different 
q vector. For monopole excitations one should use small 
wave vectors or forward scattering although a monopole 
transition contributes only to thc zero energy loss peak 
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FIG. 2: (color online) Top panels; NIXS spectra for differ- 
ent values of q in the [111] direction calculated for a NiOe 
and CoOe cluster. The monopole scattering has not been in- 
cluded. Bottom panel; NIXS intensity for the NiO peak at 3.0 
eV loss calculated for different values of q. The photon energy 
at 0=0 and the scattering angle {(p) at a photon energy of 8 
keV are given as alternative scales. 



because the excited states are orthogonal to the ground 
state since they are eigenfunctions of the same Hamiho- 
nian. Quadrupolc excitations become maximal around 
5 and hexadccapole excitations become largest be- 
tween 8 and 14 A~^. This is great as it allows one to 
choose which excitation one wants to measure. 

In the top panels of figure 2 we show the NIXS spectra 
of NiO and CoO calculated for a TMOg"" cluster consist- 
ing of a transition-metal ion and surrounding six oxygen 
ions, with the use of the program XTLS8.3 [21, [13]. For 
both NiO and CoO we see two peaks, with maximum 
intensity for q=7A~^ in good agrement with the mea- 
surements of Larson et al. The d-d excitations are 
labeled by the symmetry of the final-state without the 
inclusion of spin-orbit coupling. These peaks are split 
by spin-orbit coupling, as states of T symmetry are 3- 
fold orbital degenerate. Spin-orbit coupling has been in- 
cluded for the calculations, but the splitting can not be 
resolved with this resolution. Better experimental res- 
olutions is possible and it would be interesting to look 
at these excitations with higher resolution. We now can 
compare these calculations to the d-d spectra found in 
optical spectroscopy [l^ or RIXS 17 1 at the L2,3 edge. 
The first thing one notices is that with optics and RIXS 



one sees many more d-d excitations than with NIXS. The 
explanation is straightforward if one considers the selec- 
tion rules. For NIXS one has pure charge excitations and 
therefore the selection rule AS = 0. For NiO which has 
a ground-state of symmetry with ^296^ configuration 
[25j there are three possible excited states that are also 
triplets, namely two states of ^Ti symmetry (around 1.8 
and 3.0 eV) and one of '^T2 symmetry (around 1.1 eV). 
One can see however only two peaks as the state can 
not be reached with a quadrupolc excitation. These se- 
lection rules are rather different in RIXS at the ^2,3 edge. 
There one has an intermediate state with a core hole in 
the 2p shell of the transition metal. The spin-orbit cou- 
pling constant for 2p-core electrons of Ni is around 11.5 
eV and mixes states of different spin. This mixing re- 
sults in different spin state transitions to be observed 
with comparable intensities in RIXS at the 2p edge. 

In the bottom panel of figure 2 we show the NIXS in- 
tensity of the 3.0 eV loss peak as a function of the mag- 
nitude of q. There are several ways in which one could 
change the magnitude of q. For the geometry as shown 
in the inset of figure 2, j^fl = 2 cos(0)2i^. Which means 
one can change the energy of the photons or change the 
scattering angle in order to change the magnitude of q. 
For convenience we show three different, equivalent scales 
for the bottom panel of figure 2. 

Another advantage of NIXS is that one can not only 
tune the magnitude of q in order to optimize the scattered 
intensity one can also use the directional dependence, i.e. 
the direction of the q vector with respect to the crystal 
axes in order to do something equivalent to polarization 
analyzes. The transition matrix depends on the direc- 



tion of the q vector by, ^„ 
For a dipolc transition (fc = 1) for example this is equiv- 
alent to a dipolc in the direction of q. This allows for 
a determination of the symmetry of the initial and the 
final state, based on selection rules. In the top panels of 
figure 3 we show the angular dependence for the NIXS 
intensity of different energy loss peaks of NiO and CoO 
at 5=3.5 A^^ and q=7 . One can see that the two 
peaks of '^Ti symmetry in NiO show the same angular 
dependence whereas the two peaks in CoO, which arc of 
different symmetry show a different angular dependence. 
It should be noted that the angular dependence calcu- 
lated at (7=3.5 for the peaks at 3.0 (2.36) eV energy 
loss of NiO (CoO) show good agrement with the inten- 
sities as measured by Larson et al. [l[. It is interesting 
to note that the d-d excitation at 1.1 eV in NiO, which 
is not quadrupolc allowed can be seen at q=7 Kr^, with 
the use of a hexadecapole transition. These are strongly 
peaked in approximately the [113] direction. 

To conclude we have expanded the non-resonant con- 
tribution [A^) to the dynamical structure factor {S{q, lo)) 
in spherical harmonics. S{q, to) for CoO and NiO has 
been calculated with the use of this expansion. We used a 
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FIG. 3: (color online) Top panels; Angular dependence of the 
NIXS intensity for NiO and CoO at different loss energies 
calculated at g=3.5 and at g=7 A^^. Bottom panels; 
NIXS spectra at q=7 for different sample orientations. 



configuration interaction cluster calculation for a NiOg''" 
and CoOg°~ cluster, in order to describe the final-state 
excitons correctly. The calculated spectra are in excellent 
agreement with measurements of Larson et al. [l[ . The 
spectral representation of S{q^ co) presented here gives a 
straightforward explanation of the measured energy loss 
intensity. A big advantage is that S{q, lo) in the multipole 
expansion is easy to calculate. This is especially suitable 
for q values comparable to atomic dimensions. For larger 
energy transfers involving interband transitions and col- 
lective modes one could rely on LDA-I-U or time depen- 
dent DFT [2^, [2^. By changing the magnitude of q, one 
can tune the sensitivity of the measurement to different 
multipoles and optimize the intensity of the d-d excita- 
tion. A certain multipole has optimal intensity if the 
spherical Bessel function of the same order has a period 
comparable to the size of the local d orbital. Rotating 
the sample with respect to the q vector allows one to do 
something equivalent to polarization analysis in normal 
spectroscopy. This creates the opportunity to determine 
the symmetry of the ground-state and excited-state with 
the use of selection rules. It is important to note that 
this kind of measurement is bulk sensitive and can be 



used to study buried interfaces. The elemental sensitiv- 
ity is not as strong as in RIXS but because the radial 
matrix elements depend strongly on the radial extend of 
the d wave functions some degree of elemental sensitivity 
remains. We believe that this kind of measurement can 
provide important information on the electronic struc- 
ture and local symmetry of some of the most fascinated 
strongly correlated electron systems. 
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agencies NSERC, CIAR, and CFI. 



[1 
[2: 

[3: 

[4: 

[s: 
[6: 

[7 

[8 
[9 

[10: 
[11 

[12 

[13 
[14 

[15 

[16 

[17 

[18 
[19 



B. C. Larson, J. Z. Tischler, W. Ku, C.-C. Lee, O. D. Re- 
strepo, A. G. Eguiluz, P. Zschack, and K. D. Finkelstein, 
Subbmitted to Phys. Rev. Lett. (2007). 

N. Tsuda, K. Nasu, F. Atsushi, and K. Siratori, Elec- 
tronic Conduction in Oxides (Springer- Verlag, Berlin, 
2000). 

M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998).' 

See for examples the contributions to the Focus on Or- 
bital Physics, New Journal of Physics 6 (2004). 
Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 
J.-H. Park, L. H. Tjeng, A. Tanaka, J. W. Allen, C. T. 
Chen, P. Metcalf, J. M. Honig, F. M. F. de Groot, and 

G. A. Sawatzky, Phys. Rev. B 61, 11506 (2000). 

M. W. Haverkort, Z. Hu, A. Tanaka, W. Reicheh, S. V. 
Streltsov, M. A. Korotin, V. I. Anisimov, H. H. Hsieh, 

H. -J. Lin, C. T. Chen, et al., Phys. Rev. Lett. 95, 196404 
(2005). 

A. Ohtomo and H. Y. Hwang, Nature 427, 423 (2004). 
S. Okamoto and A. J. Millis, Nature 428, 630 (2004). 
R. Newman and R. M. Chrenko, Phys. Rev. 114, 1507 
(1959). 

E. Saitoh, S. Okamoto, K. T. Takahashi, K. Tobe, 
K. Yamamoto, T. Kimura, S. Ishihara, S. Maekawa, and 
Y. Tokura, Nature 410, 180 (2001). 

M. Griininger, R. Riickamp, M. Windt, P. Rentier, C. Zo- 
bel, T. Lorenz, A. Freimuth, and A. Revcolevschi, Nature 
418, 39 (2002). 

A. Kotani and S. Shin, Rev. Mod. Phys. 73, 203 (2001). 
J. van der Brink and M. van Veenendaal, J. Phys. Chem. 
Solids 66, 2145 (2005). 

P. M. Platzman and E. D. Isaacs, Phys. Rev. B 57, 11107 
(1998). 

C. -C. Kao, W. A. L. Caliebe, J. B. Hastings, and J.-M. 
Gillet, Phys. Rev. B 54, 16361 (1996). 

G. Ghiringhelh, M. Matsubara, C. Dallera, F. Fra- 

cassi, R. Gusmeroli, A. Piazzalunga, A. Tagliaferri, 

N. B. Brookes, A. Kotani, and L. Braicovich, Journal 

of Physics: Condensed Matter 17, 5397 (2005). 

M. Magnuson, S. M. Butorin, J.-H. Guo, and J. Nord- 

gren, Phys. Rev. B 65, 205106 (2002). 

P. Kuiper, J.-H. Guo, C. Sathe, L.-C. Duda, J. Nord- 



5 



gren, J. J. M. Pothuizen, F. M. F. de Groot, and G. A. 

Sawatzky, Phys. Rev. Lett. 80, 5204 (1998). 
[20] W. Ku, H. Rosner, W. E. Pickett, and R. T. Scalettar, 

Phys. Rev. Lett. 89, 167204 (2002). 
[21] G. J. M. Janssen and W. C. Nieuwpoort, International 

Journal of Quantum Chemistry 34, 679 (1988). 
[22] R. D. Cowan, The theory of atomic structure and spectra 

(University of Cahfornia press, Berkeley, 1981). 
[23] A. Tanaka T. Jo, J. Phys. Soc. Jpn.63, 2788 (1994). 



[24] The parameters used for the cluster calculations are 
equivalent to the ones used for the RIXS calculations 

[25] C. J. Ballhausen, Ligand field theory (McGraw-Hill book 

company. Inc., New York, 1962). 
[26] A.G. Eguiluz, O.D. Restrepo, B.C. Larson, J.Z. Tischler, 

P. Zschack and G.E. Jellison, J. Phys. Chem. Solids 66, 

2281 (2005). 



